ANALYSING POLYNUCLEOTIDE SEQUENCES 

This is a CIP of Serial No. 07/573,317 filed 
t September 28, 1990. 

1. INTRODUCTION 

Three methods dominate molecular analysis of 
nucleic acid sequences: gel electrophoresis of 
restriction fragments, molecular hybridisation, and the 
rapid ON A sequencing methods. These three methods have 
5 a very wide range of applications in biology, both in 
basic studies, and in the applied areas of the subject 
such as medicine and agriculture. Some idea of the 
scale on which the methods are now used is given by the 
rate of accumulation of DNA sequences," which is now 

10 well over one million base pairs a year. However, 

powerful as they are, they have their limitations. The 
restriction fragment and hybridisation methods give a 
coarse analysis of an extensive region, but are rapid; 
sequence analysis gives the ultimate resolution, but it 

15 is slow, analysing only a short stretch at a time. ... 
There is a need for methods which are faster than the 
present methods, and in particular for methods which 
cover a large amount of sequence in each analysis. 

This invention provides a new approach which 

20 produces both a fingerprint and a partial or complete 
sequence in a single analysis, and may be used directly 
with complex DKAs and populations of RNA without the 
need for cloning. 

In one aspect the invention provides 

25 apparatus for analysing a polynucleotide sequence, 

comprising a support and attached to a surface therof 
an array of the whole or a chosen part of a complete 
set of oligo nucleotides of chosen lengths, the 
different oligonucleotides occupying separate cells of 

30 the array and being capable of taking- part in 

hybridisation reactions. For studying differences 
between polynucleotide sequences, the invention- 
provides in another aspect apparatus comprising a 
support and attached to a surface therof an array of 

35 the whole or a chosen part of a complete set of 
oligonucleotides of chosen lengths comprising the 
polynucleotide sequences, the different 
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oligonucleotides occupying separate cells of the array 
and beingi capable of taking part in hybridisation 
reactions. 

In another aspect, the invention provides a 
5 method of analysing a polynucleotide sequence, by the 
use of a support to the surface of which is attached an 
array of the whole or a chosen part of a complete set 
of oligo nucleotides of chosen lengths, the different 
oligonucleotides occupying separate cells of the array, 

TO which method comprises labelling the polynucleotide 
sequence or fragments thereof to form labelled 
material, applying the labelled material under 
hybridisation conditions to the array, and observ i ng 
the location of the label on the surface associated 

15 with particular members of the set of oligonucleotides. 

The idea of the invention is thus to provide 
a structured array of the whole or a chosen part of a 
complete set of oligonucleotides of one or several 
chosen lengths. The array,. which may be laid out on a 

20 supporting film or glass plate, forms the target for a 
hybridisation reaction. The chosen conditions of 
hybridisation and the length of the oligonucleotides 
must at all events be sufficient for the available 
equipment to be able to discriminate between exactly 

25 matched and mismatched oligonucleotides. In the 
hybridisation reaction, the array is explored by a 
labelled probe, which may comprise oligomers of the • 
chosen length or longer polynucleotide sequences- or 
fragments, and whose nature depends on the particular 

30 application. For example, the probe nvay comprise 

labelled sequences amplified from genomic DMA by the 
polymerase chain reaction, or a mRNA population, or a 
complete set of oligonucleotides from a complex 
sequence such as an entire genome. The end result is a 

35 set of filled cells correspond i ng to the 

oligonucleotides present in the analysed sequence, and 
a set of "empty" sites corresponding to the-, sequences 
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which are absent in the analysed sequence. The pattern 
produces a fingerprint representing all of the sequence 
analysed. In addition, it is possible to assemble most 
or ail of the sequence analysed if an oligonucleotide 
5 length is chosen such that most or all oligonucleotide 
sequences occur only once. 

The number, the length and the sequences of the 
oligonucleotides present in the array "lookup table" 
also depend on the application. The array may include 
10 ail possible oligonucleotides of the chosen length, as 
would be required if there was no sequence information 
on the sequence to be analysed. In this case, the 
-referred length of oligonucleotide used depends on the 
length of the sequence to be analysed, and is such that 
15 there is likely to be only one copy of any particular 
oligomer in the sequence to be analysed. Such arrays 
are large. If there is any information available on 
the sequence to be analysed, the array may be a 
selected subset. For the analysis of a sequence which 
20 is known, the size of the array is of the same order as 
length of the sequence, and for -any applications, such 
as the analysis of a gene for mutations, it can be 
:uite small. These factors are discussed in detail in 
what follows. 
25 2. OLIGONUCLEOTIDES AS SEQUENCE ??.03ES 

Oligonucleotides form base paired duplexes with 
z 1 igonuc leo t ides which have the complementary base 
sequence. The stability of the duplex is dependent cn 
ohe length of the oligonucleotides and on base 
<0 composition. Effects of base composition on duplex 
stability can be greatly reduce: by the presence cf 
high concentrations of quaternary or tertiary amines. 
However, there is a strong effe:: of mismatches in the 
oligonucleotides duplex on tne thermal stability cf the 
' 33 r.vbrid, and it is this which ma>es the technique of 



hybridisation with oligonucleotides such a powerful 
method for the analysis of mutations, and for the 
selection of specific sequences for amplification by 
DNA polymerase chain reaction. The position of the 
mismatch affects the degree of des tab i 1 i sa tion . 

Mismatches in the centre of the duplex may cause a 

o o 
lowering of the Tm by 10 C compared with 1 C for a 

terminal mismatch. There is then a range of 

discriminating power depending on the position of 

mismatch, which has implications for the method 

described here. There are ways of improving the 

discriminating power, for example by carrying out 

hybridisation close to the Tra of the duplex to reduce 

the rate of formation of mismatched duplexes, and by 

increasing the length of oligonucleotide beyond what is 

required for unique representation. A way of doing 

this systematically is discussed. 

3. ANALYSIS OF A PREDETERMINED SEQUENCE 

One of the most powerful uses of oligonucleotide 

probes has been in the detection of single base changes 

in human genes. The first example was the detection of 

the single base change in the betagiocin gene which 

leads to sickle ceil disease. There is a need to 

extend this approach to genes in which there may be a 

number z f different mutations leading to the same 

phe no type , for example the DMD gene and the HPRT gene, 

and to find' an efficient way of scanning the human 

genome for mutations in regions which have been shown 

by linkage analysis to contain a disease' locus for 

example Huntington's disease and Cysti: Fibrosis. Any 

known sequence can be represented completely as a set 

cf overlapping oligonucleotides. The size of the set 

is li s * l ^ N, where N is the iengtr. cf the sequence 

and s is the length cf an oligomer. A gene of s kb for 

example, may be divided into an overlapping set of 
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around one tr.cusand oligonucleotides of any ca.-.en 
'■ength. An array constructed with each of cr.a» 
oligonucleotides in a separate cell can be use, as a 
multiple hybridisation probe to examine the homologous 
sequence in 'any context, a single-copy gene in the nu,an 
S enome or a messenger RNA among a mixed RNA population, 
for example. The length s may be chosen such that 
fcror - is only a small probability that any oligomer in 

sequence is represented elsewhere in the sequence 
to be analysed. This can be estimated from tne 
exoression given in the section discussing statistics 
below. For a less complete analysis it would ce 
oossible to reduce the size of the array e.g. Dy * 
factor of up to 5 by representing the sequence in a 
Partly or non-overlapping set. The advantage of using 
a comoletely overlapping set is that it provides a more 
precise location of any sequence difference, as the 
mi smatch will scan in s consecutive oligonucleotides. 
4. ANALYSIS OF AN UNDETERMINED SEQUENCE 

The genomes of all free living organisms are 
larger than a million base pairs and none has yet been 
sequenced completely. Restriction site mapping reveals 
only a small part of the sequence, and can detect only 
a small oortitn of mutations when used to compare two 
genomes. More efficient methods for analysing complex 
sequences. are needed to bring the full power =: 
molecular' genetics to bear on the many biciogtcal 
Problems for vnich there is no direct .:«« '.= the 

-r^roTvp- 4 T n many cases.* the fu.l 
gene cr genes :r.volve.. J 

< . a „ « a <- r^ei net be i e - e r ~ i n e c ; 
50 seouer.ee of tr.e nucle.c a---s 

_..,.„,.. s are those whicn differ between 

the important se-jenc-a - 

. , t.- ,{,,<, -h'-ee examoies: tr.e DNA 
two nucleic a:_ds . iQ S 1 /e 

secuences whicn are different between a wild type 
organism and :ne which carries a mutant car. lead the 
35 way to isolation of the relevant gene; similarly, the 
sequence differences between a cancer cell and its 
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normal counterpart car. reveal the cause cf transfor- 
mation; and the ?. N A sequences which differ between 
ceil types point to the functions which distinguish 
them. These problems can be opened to molecular 
5 analysis by a method which identifies sequen'ce 

differences. Using the approach outlined here, such 
differences can be revealed by hybridising the two 
nucleic acids, for example the genomic DMA of the two 
genotypes, or the m n N k populations of two cell types :c 

10 an array of oligonucleotides which represent all 

possible sequences. Positions in the array which are 
occupied by one sequence but not by the other show 
differences in two sequences. This gives the sequence 
information needed to synthesise probes which can then 

15 be used to isolate clones of the sequence involved. 

4.1 ASSEMBLING. THE SEQUENCE INFORMATION 
Sequences can be reconstructed by examining the 

result of hybridisation to an array. Any oligonucleo- 
tide of length s from within a long sequence, overlaps 

20 with two others over a length s-1. Starting from each 
positive oligonucleotide , the array may be examined for 
the four oligonucleotides to the left and the four to 
the right that can overlap with a one base displacement. 
If only one of these four oligonucleotides is found to 

25 be positive to the rigr.t, then the overlap and the ■ 
additional. base to the right determine s bases in the 
unknown sequence. The process is repeated in both 
directions, seeking unique matches with other positive 
oligonucleotides in tr.e array. Each unique match iziz 

30 £ base to the reconstructed sequence. 

4.2 SOME STATISTICS 

Any sequence cf length N can be broken down to a 
set o f ~ N overlapping sequences s base pairs in 
- e n g t h . (For double stranded nucleic acids, the 
35 sequence complexity : f = sequence of N base pairs is 
2N , because the two strands have different sequences. 
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but for the present purpose, this factor cf two is not 
significant). For oligonucleotides of length s, there 
are 4 s different sequence combinations. How big should 
s be to ensure that most oligonucleotides will be 
~ represented only once in the sequence to oe analysed, 
of complexity N? For a ran don sequence the expected 
number of s-tners which will be present in more than one 
copy is 

P>i *4'(1 — c~ x ( \ 4- X)) 

' 0 where 

A^OV-j-f l)/4' 

For practical reasons it is also useful to know 

how many sequences are related to any given s-mer by a 

• 5 single base change. Each position can be substituted 

by one of three bases, there are therefore 3s sequences 

related to an individual s-:er by a single base change, 

and the probability that any s-mer in a sequence of N 

bases is related to any other s-mer in that sequence 
— « s 

2 2 allowing one substitution is 3s x NM . The relative 
signals of matched and risza::n ed sequences will then 
depend on hew good the hybridisation conditions are in 
distinguishing a perfect mater, from one which differ by 
a single base. (If 4~ is an order of magnitude greater 

25. than U , there should only be a few, 3s/10, related to 
any oligonucleotide by one case change.} The 
indications are that the yield cf hybrid from the mis- 
matched sequence is a fraction of that formed by the 
perfect duplex. 

32 For what follows, it is assumed that conditions 

can be found which allow cligc nucleotides which have 
complements in the probe to z e distinguisr.ee from those 
w h i c h do not. 

4.3 ARRAY FORMAT, CONSTRUCT! jN AND SIZE 
35 To form an idea of the scale cf tne arrays needed 
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to analyse sequences of different complexity it is 
convenient to think of the array as = square matrix. 
All sequences of a given len S th can be represented just 
once in a matrix constructed by drawing four rows 
5 representing the four bases, followed by four similar 

columns. This produces a 4 x 4 matrix in which each of 
the 16 squares represents one of the 16 doublets. Four 
similar matrices, but one quarter the size, are then 
drawn within each of the original squares. This 
10 produces a 16 x 16 matrix containing all 256 tetra- 

nucleotide sequences. Repeating this process produces 
a matrix of any chosen depth, s. with a number of cells 
equal to 4 s . As discussed above, the cnoice of s is of 
great importance, as it determines the complexity of 
1-5 the sequence representation. As discussed below, a 
also determines the size of the matrix constructed, 
which must be very big for complex genomes. Finally, 
the length of the oligonucleotides determines the 
. hybridisation conditions and their discriminating power 
20 as hybridisation probes. 

Side of Matrix Number of 
( p i x e 1 = 1 : £■ V ~ ) Sheets of fii 
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3 65536 * x 

9 262144 
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The table shows the expected scale cf the arrays 
needed to perform the first analysis of a few genomes. 
The examples were chosen because they are genomes which 
have either teen sequenced by conventional procedures - 
the cosmic scale - , are in the process of being 
sequenced - the £. coli scale or for .which there has 

considerable discussion of the magnitude of the problem 
- the human scale. the table shows that the expected 
scale of the matrix approach is only a small fraction 
of the conventional approach. This is readily seen in 
the area cf X-ray film that would be consumed. It is 
also evident that the time taken for the analysis would 
be only a small fraction of that needed for gel 
methods. The "Genomes" column shows the length of 
random sequence which would fill about 5% of cells in 
the matrix. This has been determined to be the optimum 
condition for the first step in the sequencing strategy 
discussed below. At this size, a high proportion of 
the positive signals would represent single occurrences 
of each oligomer, the conditions needed tc compare two 
genomes for sequence differences. 
5. REFINEMENT Or AN INCOMPLETE SEQUENCE 

Reconstruction of a complex sequence produces a 
result in which the reconstructed sequence is ■ 
interrupted at any point where an oligomer that is 
repeated i r. the sequence occurs. Some repeats are 
present as components of long repeating structures 
which form part cf the structural organisation of the 
DMA , dispersed and tandum repeats in human Z N A for 
example. Eut when the length of oligonuclectide used 
in the matrix is smaller than that nee dec* tc give 
totally unique sequence representation, repeats occur 
by chance. Such repeats are likely to be isolated. 
That is, tr.s sequences surrounding the repeated 
oligomers are unrelated to each other. The gaps caused 
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by these repeats can be removed by extending the 
sequence to longer oligomers. In principle, those 
sequences shown to be repeated by the first analysis, 
using an array representation of all possible 
5 oligomers, could be r esyn the s ised with an extension at 
each end. For each repeated oligomer, there would be $ 
x 4 z 16 oligomers in the new matrix. The 
hybridisation analysis would now be repeated until tne 
sequence was complete. In practice, because the 
10 results of a positive signal in the hybridisation may 

be ambiguous, it may be better to adopt a refinement of 
the first result by extending all sequences which zLz 
not give a clear negative result in the first analysis. 
An advantage of this approach is that extending ir.e 
15 sequence brings mismatches whicn are close zo tr.e er.cs 
in the shorter oligomer, closer to the centre in tr.e 
extended oligomer, increasing the discriminatory cover 
of duplex formation. 

5.1 A HYPOTHETICAL ANALYSIS Or THE SEQUENCE OF 
20 BACTERIOPHAGE \ DNA 

Lambda phage SNA is 43,502 base pairs long. lis 
sequence has been completely determined, we have 
treated one strand :f this as a test case in a ccmputer 
simulation of the analysis. The table shows that tr.e 

25 appropriate size of oligomer to use for a sequence : f 
this complexity is the 10-mer. With a matrix of ; I- 
mers, the size was r 0 2 4 lines square. After 
"hybridisation" of the lambda 1 j-mer s in the computer , 
^ 5 ,277 ceils were positive , 1 957 had double 

30 occurrences, 75 triple occurrences, and three quadruple 
occurrences. These -5,577 positive cells representee 
known sequences, de-ermined from their position in tr.e 
matrix. Each was extended by four x one base at tr.e 3 ' 
end and four x one base at the 5 ' , end to give \i x 

35 ^6,377 = 7*12,032 cells. This extended set reduce! tr.e 



number of double occurrences to 161, a further 16- 
fold extension brought the number down to 10, and one 
more provided a completely overlapped result. Of 
course, the same end result of a fully overlapped 
sequence could be achieved starting with a 4 1 ^ matrix, 
but the matrix would.be UCOO times bigger than the 
matrix needed to represent all 10-mers, and most of the 
sequence represented on it would be redundant. 
5.2 LAYING DOWN THE MATRIX 

The method described here envisages that the 
matrix will be produced by synthesising oligonucleotides 
in the cells of an array by laying down the precursors 
for the four bases in a predetermined pattern, an 
example of which is described above. Automatic 
equipment for applying the precursors has yet to be 
developed, but there are.covious poss' ib il i t es ; it 
should not be difficult tc adapt a pen plotter or other 
computer-controlled printing device to the purpose. 
The smaller the pixel size of the array the better, as 
complex genomes need very large numbers of cells. 
However, there are limits tc how small :'nese can be 
made. 100 microns would "re a fairly comfortable upper 
limit, but could probably r. z t be achieved on paper for 
reasons of texture and diffusion. On a smooth 
impermeable surface, sucn as glass, it may be possible 
to achieve a resolution z f around 10 microns, for 
example by using a laser typesetter to preform a 
solvent repeiiant grid, building the oligonucleo- 

tides in the exposed re:i::.3 . One attractive 
possibility, which allows adaptation of present 
techniques of oligcnucler tide synthesis, is tc sinter 
microporo us glass in micrcsccpic patches onto tne 
surface of a glass plate. Laying down very large 
number cf lines or acts c::l: take a long time, if tne 
printing mechanism were slew. However, a low cost ink- 
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jet printer can print at speeds of about 10,000 spots 
per second. With this sort of speed, 10 spots could 
be printed in about three hours . 



5 5.3 OLIGONUCLEOTIDE SYNTHESIS 

There are several methods of synthesising 
oligonucleotides. Most methods in current use attach 
the nucleotides to a solid support of controlled pore 
size glass (CPG) and, are suitable for adaptation to 

10 synthesis on a glass surface. Although we know of no 
description of the direct use of oligonucleotides as 
hybridisation probes while still attached to the matrix 
on which they were synthesised, there are reports of 
the use of oligonucleotides as hybridisation probes on 

15 solid supports to which they were attached after 

synthesis. PCT Application WO 85/01051 describes a 
jc^thod for synthesising oligonucleotides tethered to a 
CPG column. In an experiment performed by us, CPG was 
used as the support in an Applied Bio-sytems 

20 oligonucleotide synthesiser to synthesise a 13-mer 
complementary to the left hand cos site of phage 
lambda- The coupling steps were all close to 
theoretical yield. The first base was stably attached 
to the support medium through all the synthesis and 

25 deprotection steps by a covalent link. 

5.4 ANALYSING SEVERAL SEQUENCES SIMULTANEOUSLY 

The method of this invention can be used to 
analyse several polynucleotide sequences 
simultaneously. To achieve this, the oligonucleotides 
may be attached to the support in the form of (for 
example) horizontal stripes. A technique for doing 
this is described in Example 3 below. Each DNA sample 
to be analysed is labelled and applied to the surface 
carrying the oligonucleotides in the form of a stripe 
(e.g. vertical) orthogonal to the oligonucleotide 
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stripes of the array. Hybridisation is seen at the 
intersections between oligonucleotide stripes and 

between them. 

5 Where sequence variations are known, an 

advantage of using this technique is that many 
. different mutations can be probed simultaneously by 
laying down stripes corresponding to each allelic 
variant. With a density of one oligonucleotide per 
10 mm, and one "individual" per 5 mm, it should be 

possible to analyse 2000 loci on a plate 'l 00 ram square. 
Such a high density of information, where the 
oligonucleotides do identify specific alleles;, is not 
available by other techniques . 

15 

6- PROBES, HYBRIDISATION AND DETECTION 

The yield of oligonucleotides synthesised on 

microporous glass is about 30 iimol/g. A patch of this 

material 1 micron thick by 10 microns square would hold 
-12 

2 q * 3 x 10 pjtiol, equivalent to about 2 g of human 
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DNA. The hybridisation reaction could therefore be 
carried out with a* very, large excess cf tr.e bound 
oligonucleotides over that in the probe. So it should 
be possible to design a system capable of 
5 distinguishing between hybridisation involving single 
and multiple occurrances of the probe sequence, as 
yield will be proportional to concentration at ail 
stages in ■ tb.^reaction . 

The polynucleotide sequence to be analysed may be 

10 of DNA or RNA. To prepare the probe, the polynucleotide 
may be degraded to form fragments. Preferably it is 
degraded by a method which is as random =3 possible, to 
an average length around the chosen ler.gtr. s of the 
oligonucleotides on the support, and oligomers of e x a : t 

15 length s selected by electrophoresis on a sequencing 

gel. The probe is then labelled. For example, 

oligonucleotides of length s may be end labelled. If 
3 2 

labelled with J ?, the radioactive yield cf any 
individual s-mer even from total human DNA could be 

20 more than 10 dpm/mg of total DMA. For detection, only 
a small fraction of this is needed in a :=:ch 10-100 
microns square. This allows hy br id i s a t i z r. conditions 
to be chosen to be close to the Tm of duplexes , which 
decreases the yield of hybrid and decreases the rate cf 

25 formation, but increases the discriminating power. 

Since the bound oligonucleotide is in excess, signal 
need not be a problem even working close to 
equilibrium. 

Hybridisation conditions can be chosen to be those 
30 known to be suitable in standard procedures used to 
hybridise to filters, but establishing optimum 
conditions is important. In particular, temperature 
needs to be controlled closely, preferably to better 
than +G.5"C. Particularly when the cncssr. length of 
25 the oligonucleotide is small, the analysis needs to be 
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able to distinguish between slight differences of rate 
and/or extent of hybridisation. The equipment may need 
to be programmed for differences in base composition 
between different oligonucleotides. In constructing 
5 the array, it may be preferable to partition this into 
sub-matrices with similar base compositions. This may 
make it easier to define the Tm which may differ 
slightly according to the base composition. 
The choice of hybridisation solvent is 

10 significant. When 1M NaCl is used, G:C base pairs are 
more stable than A:T base pairs. Double stranded 
oligonucleotides with a high G+Ccontent have a higher 
Tm than corresponding oligonucleotides with a high A+T 
content. This discrepancy can be compensated in 
various ways: the amount of oligonucleotide .laid down 
on the surface of the support can be varied depending 
on its nucleotide composition; or the computer used to 
analyse the data can be programmed to compensate for 
variations in nucleotide composition. A preferred 

20 method, which can be used either instead of or in 
addition to those already mentioned, is to use a 
chaotropic hybridisation solvent, for example; a 
quarteraary or tertiary amine as mentioned above. 
Tetramethylammoniumchloride (THAC1) has proved 

25 particularly suitable, at concentrations in the range 2 
M to 5.5 M. At TMAC1 concentrations around 3,. 5 M to 4 
M, the T^ dependence on nucleotide composition is 
greatly reduced. 

The nature of the hybridisation salt used also has 

30 a major effect on the overall hybridisation yield. 

. Thus, the use of TMAC1 at concentrations up to 5 M can 
increase the overall hybridisation yield by a factor of 
30 or more (the exact figure depending to some extent 
on nucleotide composition) in comparison with 

35 hybridisation using 1M NaCl. Manifestly, this has 
important implications; for example the amount of 
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probe material that needs to be used to achieve a given 
signal can be much lower, 

32 

Autoradiography , especially with P causes image 
degradation which may be a limiting factor 'determining 
5 resolution; 'the limit for silver halide films is around 
25 microns* Obviously some direct detection system 
would be better. Fluorescent probes are envisaged; 
given the high concentration of the target 
oligonucleotides, the low sensitivity of fluorescence 

10 may not be a problem. 

We have" considerable experience of scanning auto- 
radiographic images with a digitising scanner. Our 
present design is capable of resolution down to 25 
microns, which could readily be extended down to less 

^ than present application, depending on the quality of 
the hybridisation reaction, and how good it is at 
distinguishing absence of a sequence from the presence 
of one or more. Devices for measuring astronomical 
plates have an accuracy around 1 p. Scan speeds are 

20 such that a matrix of several million cells can be 

scanned in a few minutes. Software, for .the analysis of 
the data is straight-forward, though the large data 
sets need a fast computer. 

Experiments presented below demonstrate the 

2 g feasibility of the claims. 

Commercially available microscope slides (BDH. 
Super Premium 76 x 26 x 1 mm) were used as supports. 
These were derivatised with a long aliphatic linker 
that can withstand the conditions used for the 

30 
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deprotection of the aromatic heterocyclic bases, i.e. 

30% NH at 55° for 10 hours. The linker, bearing a 

hydroxyl group which serves as a starting point for the 

subsequent oligonucleotide, is synthesised in two 

steps. The slides are first treated with a 25% 

solution of 3-glyc idcxyprcpy 1 tri ethoxys i 1 ane in xylene 

containing several drops of Hunig's base as a catalyst. 

The reaction is carried out in a staining jar, fitted 

o 

with a drying tube, for 20 hours at 90 C. The slides 
are washed with MeOH, Et^O and air dried. Then neat 
hexaethylene glycol and a trace amount of cone, 
sulphuric acid are added and the mixture kept at 30 
for 20 hours. The slides are washed with MeOH, Et 0 

O c 

air dried and stored desiccated at -20 until use. 
This preparative technique is described in British 
Patent Application 8822228.6 filed 21 September 1988. 

The oligonucleotide synthesis cycle is performed 
as follows: 

The coupling solution is made up fresh for each 
step by mixing 6 vol. of C.5M tetrazole in anhydrous 
acetonitrile with 5 vol. cf a 0.2M solution of the 
required beta-cy ancethy 1 chosphoram i d i te . Coupling tin 
is three minutes. Oxidation with a 0.1M solution 
of I in THF/pyridine/H^O yields a stable phospho- 
triester bond. Detrity lation of the 5' end with 3% 
trichloroacetic acid in dichlororaethane allows further 
extension of the oligonucleotide chain. There was nc 
capping step since the excess cf phos:horamid i tes usee 
over reactive sites on the slide was large enough tc 
drive the coupling to ccrr.: letion. Af-ter the synthes : 
is completed, the oligonucleotide is deprotected in 30 
NH for 10 hours at 55°. The chemicals used in the 
coupling step are moisture-sensitive, and this critics 
step must be performed uncer annydrous conditions in a 
sealed container, as. follows. The shape of the pater 
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to be synthesised was cut out of a sheet of silicone 
rubber (76 x 26 x 0.5 mm) which was sandwiched between 
a microscope slide, derivatised as described above, and 
a piece of teflon of the same size and thickness. To 
this was fitted a short piece cf plastic tubi'ng that 
allowed us to inject and withdraw the coupling solution 
by syringe and to flush the cavity with Argon. The 
whole assembly was held together by fold-back paper 
clips. After coupling the set-up was disassembled and 
the slide put through the subsequent chemical reactions 
(oxidation with iodine, and det r i ty 1 at i on by treatment 
with TCA) by dipping it into staining jars. 
EXAMPLE 1. 

As a first example we synthesised the sequences 
cl igo-dT -ol igo-dT on a slide by gradually decreasing 
the level of the coupling solution in steps 10 to 14. 
Thus the 10-mer was synthesised cn the upper part of 
the slide, the 14-mer at the "bottom and the 11, 12 and 
13-mers were in between. We^used 10 pmol oligo-dA 
20 labelled at the 5 1 end with P by the polynucleotide 
kinase reaction to a total activity of 1.5 million 
c.p.m., as a hybridisation probe. Hybridisation was 
carried out in a perspex (Plexiglas) container made to 
■it a microscope slide, filled with 1.2 ml cf 1M NaCl 
25 in TE, 0.1% SOS, for 5 minutes at 20 . After a..short 
rinse in the same solution with:*Jt oligonucleotide, we 
were able to detect more than 20 ~0 c.p.s. with a 
radiation monitor. An autorac 1 : graph, .showed that all 
the counts came from. the area were the oligonucleotide 
~ a d been synthesised, i.e. there was no non-specific 
binding to the glass or to the r e g i o n that had been 
cerivatised with the linker only. After partial 
elution in 0.-1 M NaCl differential binding to the 
target is detectable, i.e. less binding to the shorter 
tnan the longer oligo-dT. By gradually heating the 
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slide in the wash solution we determined the T (mid- 

om 

point o r transition when 50* eluted) tc be 33 . There 

o 

were no counts detectable after incubation at 39 . 
The hybridisation and melting was repeated eight tirr.es 
with no diminution of the signal. The result is 
reproducible. We estimate that at least 5% of the 
input counts were taken up- by the slide at each cycle. 
EXAMPLE 2. 

In order to determine whether we would be able to 
distinguish between matched and mismatched oligo- 
nucleotides we syrvthesised two sequences 3 1 CCC GCC GCT 
GGA (cosL) and 3' CCC GCC TCT GGA, which' differ by one 
base at position 7. All bases'except the seventh were 
added in a rectangular patch. At the seventh base., 
lb half of the rectangle was exposed in. turn to add the 

two different bases, in two stripes. Hybridisation of 
cosR prcbe oligonucleotide (5' GGG CGG CGA CCT) (kinase. 

labelled with P to 1.1 million c.p.m., 0.1 M NaCl , 

o 

TE, C.15 SDS) was for 5 hours at 32 . The front of the 

^° slide showed 100 c.p.s. after rinsing. Autoradiography 
showed that annealing occurred only t: "he part of the 
slide with the fully como 1 ement ary o 1 i z :nuc 1 eot i de . 
No signal was detectable on the patch with the 
mismatcred sequence. 

25 EXAMPLE 3. 

For a further study of the effects cf mismatches 
or shorter sequences on hybridisation behaviour, we 
constructed two arrays; one (a) of 2* ,c \ igonucleot ides 
and the ether (b) of 72 oligonucleotides. 

3° These arrays were set out as shew n in Table 1(a) 

and 1(b'. . The masks used to lay down t~ese arrays were 
differs r~ from those used in previous experiments. 
Lengths cf silicone rubber tubing (1mr. c.d.) were glued 
with silicone rubber cement to the surface of plain 

35 mi c rose:: a slides, in the form of a "U" . Clamping 
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these masks against a derivatised microscope slide 
produced a cavity into which the coupling solution was 
introduced through a syringe. In this way cnly the 
part of the slide within the cavity came into contact 
with the phosphorami d i te solution. Except in the 
positons of the mismatched bases, the arrays listed in 
Table 1 were laid down using a mask which covered most 
of the width of the slide. Off-setting this .task by 
3mm up or down the derivatised slide in subsequent 
coupling reactions produced the ol ignucleotides 
truncated at the 3' or 5' ends. 

For the introduction of mismatches a -mask was used 
which covered half (for array (a)) or one third (for 
array (b)) of the width of the first mask. The bases 
15 at positions six and seven were laid down in two or 

three longitudinal stripes. This led. to the synthesis 
of oligonucleotides differing by one base on each half 
(array (a)) or third (array (b)) of the slide. In 
other positions, the sequences differed from the 
longest sequence by the absence of bases at the ends. 

In array (b!, there were two columns c r sequences 
between those shown in Table 1(b), in which the sixth 
and seventh bases were missing in all positions, 
because the slide was masked in a stripe by t u e 
silicone rubber seal. Thus there were a tctal of.. 72 
different sequences represented on the slide i~ 90' 
d i f f erent positions. 

The 19-mer 5' CTC CTG AGG AGA AGT CTG C n : s used 
for hybridisation (2 million cpm, 1.2 "ml O.'V \ = C1 in 
3° TE, 0.1% SDS, 20 C ). 

The washing and elution steps were relieve:! by 
autoradiography. The slide was kept in the washing 
solution for 5 rn 1 n at each elution step and tren 
exposed (45 ir.ir. , intensified). Elution terr.pe^etures 
" were 23, 36, *2, 47, 55 and 60 C respectively. 
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As indicated in the table, the oligonucleotides 

showed different meltina behaviour. Short oligo- 

. " o 
nucleotides melted before longer ones, and at 55 C, 

only the perfectly matched 19-mer was stable, all other 

5 oligonucleotides had been eluted. Thus the method can 

differentiate between a 18-mer and a 19-mer which 

differ only by the absence of one base at the end. 

Mismatches at the end of the oligonucleotides and at 

internal sites can all be melted under conditions where 

10 the perfect duplex remains. 

Thus we are able to use very stringent 

hybridisation conditions that el iminate .-anneal ing to 

mismatch sequences or to oligonucleotides differing ir. 

length by as little as one base. No other method using 

' 5 hybridisation of oligonucleotides bound to the solid 

supports is so sensitive to the effects of mismatching, 

EXAMPLE 4. 

To test the application of the invention to 
diagnosis of inherited diseases, we hybridised the 

20 array (a), which carries the oligonucleotide sequences 

specific for the wild type and the sickle cell muta-;:ns 
of the P -globin gene, with a 110 base pair fragment : r 
ON A amplified from- the p -clobin gene by means of the 
polymerase chain reaction (PCR). Total DNA from the 

25 blood of a normal individual (1 microgram) was 

amplified by PCR in the presence of apprcpri ate primer 
oligonucleotides. The resulting 110 base pair fragment 
was purified by elecrcchcresis through an agarose gel. 
After elution, a small sample (ca. 10'piccgram) was 

30 labelled by using << -"?-dCTP (50 microCurie) in a 
second PCR reaction. This PCR contained only the 
upstream priming ol iccr.ucleotide. After 50 cycles c s 
amplification with an extension tine cf 9 min. the 
product was removed fr:r. precursors by gel filtrati". 

35 Gel electrophores i s • z £ the radioactive product showec i 



major band corresponding in length to the 110 base 

fragment. One quarter of this product (100,000 c.p.m. 

in 0.9 M NaCl, TE, 0,1% SOS) was hybridised to the 

o 

array (a). After 2 hours at 30 ca. 15000 c.p.m. had 
been taken uo. The melting behaviour of the hybrids 
was followed as described for the 19-mer in example 3, 
and it was found that the melting behaviour was similar 
to that of the oligonucleotide. That is to say, the 
mismatches considerably reduced the melting temperature 
of the hybrids, and conditions were readily found such 
that the perfectly matched duolex remained whereas the 
mismatched duplexes had fully relted. 

Thus the invention can be used to analyse long 
fragments of DNA as well oligonucleotides, and this 
example shows how it may be used to test nucleic acid 
sequences for mutations. In particular it shows how it 
may be applied to the diagnosis of genetic diseases. 

EXAMPLE 5. 

To test an automated system for laying down the 
precursors, the cosL oligonucleotide was synthesised 
with 11 of the 12 bases added in the way described 
above. For the addition of tre seventh base, however, 
the slide was transferred intc an Argon filled chamber 
containing a pen plotter. The pen of the plotter had 
been replaced by a component, fabricated from Nylon, 
which had the same shape and dimensions as the pen, but 
which carried a polytetraf lucrrethy lene (PTFE) tube, 
through which chemicals could : :e delivered to the 
surface of the glass slice wh::h ley on the bed of the 
plotter. A microcomputer was 'jsed zo control the 
plotter and the syringe pump »rich delivered the 
chemicals. The p_en, carryin; zhe delivery tube from 
the syringe, was moved into rrrition above the slide, 
t*e pen was lowered and the p-" activated to lay down 
coupling solution. ■"*■ Filling z~ e pen successively with 
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G, T and A phosphoramidite solutions an array of twelve 
spots was laid down in three groups of four, with three 
different oligonucleotide sequences. After 
, hybridisation to cosR, as described in Example 2, and 
5 ' autoradiography, signal was seen only over the four 
spots of perfectly matched oligonucleotides, where the 
dG had been added. 
EXA M P LE 6 

This example demonstrates the technique of 

10 analysing several DNA sequences simultaneously. Using 
the technique described in Example 3, a slide was 
prepared bearing six parallel rows of oligonucleotides 
running along its length. These comprised duplicate 
hexadecamer sequences corresponding to antisense 

15 sequences of the p-globin wild-type (A), sickle cell - 
(S) and C mutations. 

Clinical samples of AC, AS and SS DNA were 
procured. Three different single-stranded probes of 
110 nt length with approx. 70,000 c.p.m, in 100 jal 1M 

20 NaCI, TE pH 7.5, 0.1% SDS, viz AC, AS y and SS; DNA were 
prepared. Radiolabeled nucleotide was included in the 
standard PCR step yielding a double -stranded labelled 
fragment. It was made single-stranded with 
Bacteriophage X exohuclease that allowed to selectively 

2g digest one strand bearing a 5' phosphate. This was 

made possible by phosphorylating the downstream primer 
with T4 Polynucleotide kinase and ("cold") ATP prior to 
PCR. These three probes were applied as three stripes 
orthogonal to the surface carrying the six 

3Q oligonucleotide stripes. Incubation was at 3b*C for 2 
hours in a moist chamber. The slide was then rinsed at 
ambient temperature, then 45 - C for 5 minutes and 
exposed for 4 days with intensification. The genotype 
of each clinical sample was readily determined from the 
- autoradiographic signals at the points of intersection. 
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EXAMPLE 7 , 

A plate was prepared whose surface carried an 
array of all 256 octapurines. That is to say ;/ the 
array comprised 256 oligonucleotides each consisting 
5 of a different sequence of A and G nucleotides. This 
array was probed with a mixture comprising all 256 
octapyrimidines, each end labelled by means of 
polynucleotide kinase and )/- 32 P-ATP. Hybridisation 
was performed for 6-8 hours at 4*C. 

10 In consecutive experiments the hybridisation 

solvent was changed through the series 1M NaCl 
(containing 10mM Tris.HCl pH 7.5, 1mM EDTA, 7j% 
sarcosine) and 2M, 2.5M, 3M, 3.5M, 4M, 4.5M, 5M and 
5.5M TMAC1 (all containing 50mM Tris^HCl pH 8.0, 2mM 

15 EDTA, SDS at less than 0.04 mg/ml) • The plate was 

rinsed for 10 minutes at 4'C in the respective solvent 
to remove only loosely matched molecules, sealed in a 
plastic bag and exposed to a Phorphorlmager storage 
phosphor screen at 4'C overnight in the dark. 

20 The following table quotes relative signal 

intensities, at a given salt concentration, of hybrids 
formed with oligonucleotides of varying a content* in 
this table, the first row refers to the oligonucleotide 
GGGGGGGG, and the last row to the oligonucleotide 

25 AAAAAAAA. It can be seen that the difference in 

response of these two oligonucleotides is marked in 1M 
NaCl, but much less marked in 3M or 4M TMAC1. 
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Relative Intensi ti as at civap Salt C oncentration 



Solvent Number of A's 

0 4 8 



100 30 20 

100 70 30 

. 70 100 40 

60 100 40 



The following table indicates relative signal 
15 intensities obtained, with octamers containing 4A*s and 
4G's, at different hybridisation salt concentrations. 
It can be seen that the signal intensity is 
dramatically increased at higher concentrations of 
TMAC1. 

20 

Relative Intensities a t different Salt nonceiitrat-^^ 



1M NaCl 
2M TMAC1 
3M TMAC1 
4M TMAC1 



25 



Solvent Yield of hybrid 



1M NaCl .100 

2M TMAC1 200 

30 3M TMAC1 700 

4M TMAC1 '2000 



In conclusion, we have demonstrated the following: 
1 . It is possible to synthesise oligonucleotides in 
good yield on a flat glass plate. 



-26- 

'2. Multiple sequences can be synthesised on! the sample 
in small spots, at high density, by a simple manual 
procedure, or automatically using a computer controlled 
device. 

5 3. Hybridisation to the oligonucleotides ori the plate 
can be carried out by a very simple procedure;* 
Hybridisation is efficient, and hybrids can be detected 
by a short autoradiographic exposure ♦ 

4. Hybridisation is specific. There is no detectable 
qo signal on areas of the plate where there. are Ino 

oligonucleotides. We have tested the effects of mis- 
matched bases, and found that a single mismatched base 
at any position in oligonucleotides ranging in length 
from 12-mer to 19-mer reduces the stability o s f the 
15 hybrid sufficiently that the signal can be reduced to a 
very low level, while retaining significant 
hybridisation to the perfectly matched hybrid. 

5. The oligonucleotides are stably bound td the glass 
and plates can be used for hybridisation repeatedly. 

20 Ttie invention thus provides a novel way of 

analysing nucleotide sequences, which should find a 
wide range of application. We list a number of 
potential applications below: 

25 
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Small arrays of oligonucleotides as fingerprinting 

and mapping tools 

Analysis of known mutations including genetic diseases. 

Example 4 above shows how the invention may be 
used to analyse mutations. There are many 
applications for such a method, including the detection 
of inherited diseases. 

Genomic fingerprinting. 

In the same way as mutations which lead to disease 
can be detected, the method could be used to detect 
point mutations in any stretch of DNA. Secuences are 
now available for a number of regions containing the 
base differences which lead to restriction fragment 
length polymorphisms (RFLPs). An array of oligo- 
nucleotides representing such polymorphisms could be 
made from pairs of oligonucleotides representing the 
two allelic restriction sites. Amplification of the 
sequence containing the RFLP, followed by hybridisation 
to the plate, would show which alleles were present in 
the test genome. The number of oligonucleotides that 
could be analysed in a single analysis could be quite 
Large. Fifty pairs made from selected alleles would 
be enough tc give a fingerprint unique to an individual. 

Linkage analysis. 

Applying the method described in the last 
paragraph t: a pedigree would pinpoint reco.-.: inations . 
Each pair o r spots in the array would give the 
information that is seen in the track c~ the RFLP 
analysis, using gel electrophoresis and b 1 : 1 1 i n g , that 
is now routinely used for linkage studies. it should 
be possible to analyse many alleles in 2 sincle 
analysis, by hybridisation to an array of allelic pairs 
of oligonucleotides, greatly simplifying the methods 
used tc fin: linkage between a DNA polymorphism and 
phenotypic r.arker such as a disease gene. 

The exar.ples above could be carried •c'i , t using the 
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method we have developed and confirmed by experiments. 
Large arrays of oligonucleotides as sequence 

reading tools. 

We have shown that oligonucleotides can be 
synthesised in small patches in precisely determined 
positions by one of two methods: by delivering the 
precursors through the pen of a pen-plotter, or. by 
masking areas with silicone rubber. It is obvious how 
a pen plotter could be adapted to synthesise large 
arrays with a different sequence in each position. For 
some applications the array should be a predetermined, 
limited set; for other aoo 1 i cat ions , the array should 
comprise every sequence of a predetermined length. The 
masking method can be used for the latter by laying, 
down the precursors in a mask which .produces 
intersecting lines. There are many ways in which -his 
can be done and we give one example for illustration: 

1. The first four bases, A, C, 6, T, are laid in four 
broad stripes on a square plate. 

2. The second set is laid down in four stripes equal 
in width to the firs*, and orthogonal to them. The 
array is now composed of all sixteen d inuc 1 eot ides . 

3. The third and fourth layers are laid down in f:ur 
sets of four stripes one quarter the width of the firs- 
stripes. Each set cf four narrow stripes runs ,wi;Mn 
one of the broader stripes. The arrdy is now composed 
of all 256 tetranuclect ides . 

4. The process is. reoeated, each time laying down tw 
layers with stripes which are one quarter the wid-* of 
the previous two layers. Each layer added increases 
the length of the oligonucleotides by one base, an: th 
number of different oligonucleotide sequences by a 
factor of four. 

The dimensions of such arrays are determined by 
the width of the stripes. The narrowest .stri pe we 
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have laid is 1mm, but this is clearly not the lowest 
limit. 

There are useful applications for arrays in which 
part of the sequence is predetermined and part made up 
of all possible sequences. For example: 

Characterising mRNA populations. 

Most mRNAs in higher eukaryotes have the sequence 
AAUAAA close to the 3' end. The array used to analyse 
mRNAs would have this, sequence all over the plate. To 
analyse a mRNA population it would be hybridised to an 
array composed of all sequences of the type N AATAAAN 
For rn + n = 8, which should be enough to give"^ unique 
oligonucleotide address to most- of the several thousand 
mRNAs that is estimated to be present in a source such 
as a mammalian cell, the array would be 256 elements 
square. The 256 x 256 elements would be laid on the 
AATAAA using the masking method described above. With 
stripes of around 1mm, the array would be ca. 256mm 
square . 

This analysis would measure the complexity of the 
mRNA population and could be used as a basis for 
comparing populations from different cell types. The 
advantage of this approach is that the differences in 
the hybridisation pattern would provide the sequence of 
oligonucleotides that could be used as probes to 
isolate all the mRNAs that differed in the populations. 

Sequence determination. 

To extend the idea to determine unkown sequences, 
using an array composed of all possible oligonucleotides 
of a chosen length, requires larger arrays than we have 
synthesised to date. However, it is possible to scale 
down the size of spot and scale up the numbers to those 
required by extending the -ethods we have developed and 
tested on small arrays. Cur experience shows that the 
method is much simpler in operation than the gel based 
methods . 
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For Exa: 
20 GAG 



JO 

36 
47 
60 
56 
56 
47 
42 
36 
36 
36 



GAG 
GAu 

r a 

urvu 

GAG 
AG 



:ples 3 and 4 array ( 
GAC TCC TCT ACG 
GAC TCC TCT GAC G 
GAC TCC TCT AGA CG 
GAC TCC TCT CAG ACG 
GAC TCC TCT TCA GAC 
GAC TCC TCT TCA GAC 
GAC TCC. TCT TCA GAC 
GAC TCC TCT TCA GAC 
.AC TCC TCT TCA GAC 
..C TCC TCT TCA GAC 
... TCC TCT TCA GAC 
... .CC TCT TCA GAC 



TABLE 1 

a) was set cut as fol 



20 
20 
20 
36 
47 
42 
42 
42 
36 
36 
36 
36 



GAG GAC aCC 

GAC GAC aCC 

GAC GAC aCC 

GAG GAC aCC 

GAG GAC aCC 

AG GAC aCC 

G GAC aCC 

GAC aCC 

.AC aCC 

..C aCC 

• • • dCC 

• • • • CC 



lows : 
TCT ACG 
TCT GAC 
TCT AGA 
TCT CAG 
TCT TCA 
TCT TCA 
TCT TCA 
TCT TCA 
TCT TCA 
TCT TCA 
TCT TCA 
TCT TCA 



G 

CG 
ACG 
GAC G 
GAC G 
GAC G 
GAC G 
GAC G 
GAC G 
GAC G 
GAC G 
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20 



25 
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Fo>* ■ 
20 GAG 
20 GAG 
20 GAG 
20 G*G 
20 GAC 
20 GAG 
20 G^G 
20 GAC 
32 GAC 
22 GAC 
42 GAG 
52 GAG 
42 .AG 
42 ..G 
37 ... 
32 ... 
32 
22 

calu- 
eve-.. 
(a.t: 



xample 3 array (b) was set out as follows: 



GAt TC 
GAt TCC 
GAt TCC T 
•GAt TCC TC 
GAt TCC TCT 
GAt TCC TCT T 
GAt TCC TCT TC 
GAt TCC TCT TCA 
GAt TCC TCT TCA G 
GAt TCC TCT TCA GA 
GAt TCC TCT TCA GAC 



20 GAG GAC TC 
20 GAG GAC TCC 
20 GAG GAC TCC T 
20 GAG GAC TCC TC 
20 GAG GAC TCC TCT 
20 GAG G^C TCC TCT T 
20 GAG GAC TCC TCT TC 
20 GAG GAC TCC TCT TCA 
42 GAG GAC TCC TCT TCA G 
47 GAG GAC TCC TCT TCA 
52 GAG GAC TCC TCT TCA GAC 



20 GAG GAC aC 

20 . GAG GAC aCC 

20 GAG GAC aCC T 

20 GAG GAC aCC TC 

20 GAG GAC aCC TCT 

20 GAG GAC aCC TCT T 

20 GAG G£ aCC TCT TC 

20 GAG GAC aCC TCT TCA 

20 GAG GAC aCC TCT TCA 

22 GAG GAC aCC TCT TCA 



42 GAG QfiC aCC TCT TCA 
GAt TCC TCT TCA GAC G 60 GAG GAC TCC TCT TCA GAC G 52 GAG GAC aCC TtT TCA 
GAt TCC TCT TCA GAC G 52 .AG GAC TCC TCT TCA GAC G 42 .AG GAC aCC TCT TCA 
GAt TCC TCT TCA GAC G 52 ..G GAC TCC TCT TCA. GAC G 42 ..G GAC aCC TCT TCA 
GAt TCC TCT TCA G^C G 47 ... GAC TCC TCT TCA GAC G37 ... GAC aCC TCT TCA 
..At TCC TCT TCA GAC G-42 ; ; v .AC TCC TCT TCA GA! G 32 . . . .AC aCC TCT TCA 

..t TCC TCT TCA GAC £42 C TCC TCT TCA GAC G' 32 C aCC TCT TCA 

... TCC TCT TCA GAC G 32 TCC TCT TCA GAZ G 32 aCC TCT TCA 

~n the three columns of array (b) 1 i ".ed above, were two 
•. n .s, in which bases 6 and 7 from the left were missing in 
line. These sequences all melted at 20 or 32 degrees, 
mismatch base (.) missing base. 



GA 
GAC 

GAC G 

QAC G 

GAC G 

GAC G 

GiC G 

GAC G 

GAC G 



